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O " Abstract 

The connection of maximally supersymmetric Yang-Mills theory to the (2, 0) the- 
ory in six dimensions has raised the possibility that it might be perturbatively 
ultraviolet finite in five dimensions. We test this hypothesis by computing the 
coefficient of the first potential ultraviolet divergence of planar (large N c ) max- 
imally supersymmetric Yang-Mills theory in D = 5, which occurs at six loops. 
We show that the coefficient is nonvanishing. Furthermore, the numerical value 
of the divergence falls very close to an approximate exponential formula based on 
the coefficients of the divergences through five loops. This formula predicts the 
approximate values of the ultraviolet divergence at loop orders L > 6 in the critical 
dimension D = 4 + 6/L. To obtain the six-loop divergence we first construct the 
planar six-loop four-point amplitude integrand using generalized unitarity. The 
ultraviolet divergence follows from a set of vacuum integrals, which are obtained 
by expanding the integrand in the external momenta. The vacuum integrals are 
integrated via sector decomposition, using a modified version of the FIESTA pro- 
gram. 
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I. INTRODUCTION 

Recent years have seen impressive progress in computing scattering, amplitudes in general 
gauge and gravity theories (see, for example, the recent reviews The progress has 

been especially great for maximally supersymmetric Yang- Mills theory (MSYM), a theory 
with sixteen supercharges whose D = 4 version is M = 4 super- Yang-Mills theory. One 
application has been to study the ultraviolet (UV) properties of both gauge and gravity 
theories. The all-loop ultraviolet finiteness of M = 4 super- Yang-Mills theory in D = 4 
was established in the 1980's Q. In dimensions D > 4, explicit unitarity-based amplitude 
computations in MSYM in the 1990's @, 3 showed that its degree of convergence was a bit 
better than had been anticipated. The results suggested that the correct finiteness bound 
for MSYM in D dimensions at L loops is 

D<4+y (L>2). (1.1) 

This bound is consistent with all-loop ultraviolet finiteness in D = 4. The bound (11. ip has 
been confirmed to all loop orders [Bj using harmonic superspace The case L = 1 is 
an exception; at one loop the first divergence is in D = 8, not D — 10 {2|. Interestingly, 
maximal Af = 8 supergravity follows precisely the same finiteness bound, at least through 
four loops [8- Tot . 



An important remaining question is whether the bound is saturated or not; that is, 
whether the coefficient of the potential logarithmic divergence in the critical dimension 
-D = 4 + 6/Lis nonzero or not, for each loop order. On the one hand, if the theory contains 
some unknown or hidden symmetry, then the coefficient in D = 4 + Q/ L could vanish, leading 
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to a higher critical dimension than expected at some loop order. That is, at a loop order 
affected by the symmetry, the lowest dimension with an ultraviolet divergence would be 
surprisingly high. On the other hand, if the bound is saturated, it proves that no additional 
hidden symmetries exist that alter the degree of divergence — at least through the loop 
orders explored. 

The only known reliable means for answering such a question is to explicitly evaluate 
the ultraviolet divergence of an appropriate on-shell multi-loop scattering amplitude in the 
expected critical dimension. Such a computation can be performed either within the large- 
N c , or planar, limit of the theory with gauge group SU(N C ), or for a general gauge group 
inclu ding all sub leading terms in the 1/N C expansion. We know from these computations [4j, 
0, HHJ that the bound (11. ip is indeed saturated through at least five loops. Interestingly, 



certain subleading-in-iV c terms (ejj. double-trace terms), do have an improved behavior at 



three loops and beyond [10|, llll, m ^ ms P a P er ; we wm only consider the leading-color 
(planar) terms. 

If the bound (11.11) is saturated, then eq. ( II .ip implies that in D = 5 a divergence should 
first appear at six loops. A primary purpose of this paper is to compute the coefficient of 
this divergence. The case of D = 5 is especially interesting because this theory has a UV 
completion, the (2, 0) theory in D = 6. The (2, 0) theory has no Lagrangian description; 
rather its existence follows from arguments in string theory and M theory [15] . This connec- 
tion suggests that UV divergences in D = 5 could give us information about the (2, 0) theory 



16|. Of course, a low-energy effective theory usually has UV divergences, and by itself this 
connection does not lead to constraints. However, the present example is somewhat unique 
in that the (2, 0) theory has neither a dimensionless coupling constant nor a preferred scale, 
so that seemingly different D = 5 quantities turn out to be related nonperturbatively. For 
example, the Kaluza-Klein modes in D = 5, arising from the compactification of the (2, 0) 
theory on a circle, can also be identified with solitons in the gauge theory [l6l-18] (solutions 
correspon ding to instantons in D = 4). 

In ref. |l6| these aspects were discussed in the context of ^-duality of the D = 4 J\f = 4 
super- Yang- Mills theory obtained by compactification of D = 5 MSYM on a circle, which 
can also be thought of as compactification of the (2, 0) theory on a two-torus. In this 
construction S-duality has a geometric origin as an exchange of the two sides of the torus. 
This argument can be re-expressed in D = 5 terms, and UV divergences in D = 5 can 
potentially violate the D = 4 S'-duality. An alternative argument suggesting finiteness [l8[ 
is based on the soliton-Kaluza-Klein correspondence for the compactified (2, 0) theory in 
the phase where the gauge symmetry is broken by separating the multiple branes used in 
its construction. Although these arguments do not prove that there are no UV divergences, 
they do motivate the question. In the present paper we definitively answer the question, 
by computing the numerical coefficient of the potential divergence in planar MSYM at six 
loops. 

We find that the bound (II. ip is indeed saturated for L = 6 and D = 5; that is, the 
divergence has a nonzero coefficient. Somewhat surprisingly, we also find that, through at 
least six loops, the numerical values of the leading-color planar critical-dimension divergences 
can be fit accurately to a simple exponential Ansatz. Although we do not understand the 
origin of this simple functional form, it does have useful consequences: It gives us additional 
confidence that we have correctly computed the six-loop divergence and that it is non-zero. 
Moreover, by extrapolating to higher-loop orders it allows us to predict the approximate 
numerical values of the divergences for L > 7 in their critical dimensions. This result 
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suggests that even outside of four dimensions, MSYM has a surprisingly simple structure, 
reflected in the simple pattern of ultraviolet divergences. 

Until recently, a direct evaluation of the six-loop ultraviolet properties of MSYM would 
have been out of reach. However, a combination of advances has made it possible. Many 
breakthroughs in understanding the structure of integrands for multi-loop amplitudes now 
allow for a rather straightforward construction of the integrand for the six-loop four-point 
amplitude in planar MSYM In addition to the integrand construction presented here, 



two recent papers give the same integrands in D = 4 [19|, |20| , albeit presented somewhat 
differently. Our construction is also valid for loop momenta spanning the full five dimensions. 
Well-developed techniques for extracting UV divergences from amplitudes 

mm mil then 

allow us to express the six-loop divergence in terms of a relatively small number of vacuum 
integrals. Furthermore, a set of integral consistency relations [23], related to integration by 



parts identities [24J , allows us to reduce the number of integrals further. It also provides for 
nontrivial cross-checks on numerical evaluations of the integrals and an independent means 
for estimating numerical integration uncertainty. 

The final integrals obtained, after imposing the consistency relations, are nevertheless 
quite challenging to evaluate analytically. Instead, we make use of the major advances 
in numerical integration techniques. A long-standing challenge of numerically evaluating 
Euclidean Feynman integrals has been addressed by a computational technique called sector 
decomposition {25[, implemented in several software packages (26l-28| including the FIESTA 
package [29|, used here. To our knowledge, the sector decomposition technique has not 
previously been applied at such a high loop order. However, when applied to this problem it 
leads successfully to integrals that can be evaluated on a moderate size (1000 node) cluster 
in a few days. 

The remainder of this paper is organized as follows. In section[TT]we outline the derivation 
of the integrand. Section III II describes the procedure for extracting the UV divergence 
from the amplitude in terms of a set of vacuum integrals. Section IIVI explains the sector 
decomposition method, and section |V] describes details of the numerical evaluation. In 
section |VI] we give our conclusions, and comment on the feasibility of evaluating the seven- 
loop integrals that might be required to check whether the coefficient of the first potential 
counterterm of Af = 8 supergravity vanishes. 



II. CONSTRUCTING THE INTEGRAND 

Our study of the UV properties of the planar MSYM amplitude in D = 5 begins by 
constructing the integrand of the six-loop four-point amplitude. Because we need the inte- 
grand in five dimensions, we must ensure that our construction is valid for loop momenta 
inhabiting five spacetime dimensions. (We can always take the external momenta to be four- 
dimensional, and assign helicities to the external gluons if desired.) The five-dimensional 
validity of the integrand is accomplished by verifying unitarity cuts in higher dimensions, 
which we have done on a large class of cuts. In addition, to extract the UV divergences, we 
prefer a local form for the integrand, in which the only denominator factors are standard 
Feynman propagators. To find the desired representation, we use generalized unitarity, a 



particularly effective general purpose refinement of the unitarity method [30|. (For recent 



reviews of this method see refs. [3l[.) Our form for the integrand differs somewhat from 



recent ones based on four- dimensional constructions 19, l20J. However, we have confirmed 



analytic agreement in any dimension with the form in ref. [20| (which is also known to agree 
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with that in ref. [19|]). 

We will focus on the leading-color, planar contribution to the six-loop amplitude in 
SU (N c ) gauge theory, which has the same color structure as the tree amplitude, up to overall 
factors of the number of colors, N c . The color-decomposed form of the planar contribution 
to the L-loop four-point amplitude is, 



4 L) = 9 2 



g 2 N c j Tr{T a ^T a ^T a ^T a ^)A[ L) {p{l),p{2),p{3),p{4)) , (2.1) 



where is an L-loop color-ordered partial amplitude. The sum runs over non-cyclic 
permutations, p, of the external legs. In this expression we have suppressed momentum and 
helicity labels, leaving only the indices identifying the external legs. This decomposition 
holds for any set of external particles from the full gauge supermultiplet. 

We will not describe our construction of the six-loop amplitude in any detail; it is similar 
to the construction of the five-loop planar amplitude given in ref. [l2[. Integrands of planar 
amplitudes in MSYM are relatively simple to obtain because dual conformal symmetry 



severely restricts their form [12|, |32H35| . Although dimensional regularization breaks dual 
conformal invariance, it does so mildly at the level of the integrand. Indeed, the integrands 
of loop amplitudes in planar MSYM are known to have the same simple properties under 
dual conformal transformation in dimensions D < 6 (and likely in all dimensions D < 10) as 



they have in four dimensions [36H38| . The only breaking of dual conformal invariance comes 
from the integration measure. This allows us to use dual conformal symmetry to guide the 
construction of the integrand, even outside of four dimensions. 
We write the Ansatz for the six-loop planar amplitude as 

/ 10 C\ D T) 
ri(2^ x ' (2 - 2) 

where v44 ree (l, 2, 3, 4) is the color-ordered tree amplitude and the Mandelstam invariants are 
s = (ki + k 2 ) 2 and t = (k 2 + k 3 ) 2 . For bookkeeping purposes we organize the integrand in 
terms of graphs with only cubic vertices. We incorporate any contact (four-point) interac- 
tions by including numerator terms that can potentially cancel propagators. Thus there is 
no loss of generality in using cubic graphs. We decompose the integrand X as 

68 68 
D 4 i=l ° l D 4 i=l ° l 11^=5^ 

The sum runs over a set of distinct planar cubic graphs, which contribute in all eight possible 
arrangements generated by the dihedral group D4 (corresponding to symmetries of a square 
with corners labeled by the four external momenta). The dimension of the symmetry group 
leaving each diagram invariant is Si. At six loops, there are 68 non- vanishing topologically 
distinct graphs shown in Figs. Q] and [2j The product in the denominator of eq. (12.31) runs 
over the 19 internal Feynman propagators of each labeled graph. The numerators, iVj of 
each integral are polynomials, 

.V, Y."'J }, '>- ( 2 - 4 ) 
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FIG. 1: Graphs 1 through 35 for the planar six-loop four-point amplitude. 



where the monomials M^- depend only on Lorentz invariants constructed from the dual 
(loop) momenta for each diagram, and the a^- are numerical coefficients (rational numbers) 
to be determined from various constraints. 

As a first step, we require the monomials to have the proper weight under dual conformal 
transformations. To expose the dual conformal properties we use the standard |32| dual 
variables X{ — Xj = Xij, with 



241 — k\ , X\2 — &2> ^23 — ^3; 234 — £4 



(2.5) 



where ki are the external momenta. As discussed in detail in ref. (I2j, a practical way of 
expressing the internal momenta of a diagram in terms of dual variables is to use an (L+ 1)- 
particle cut, which divides the L-loop amplitude into two tree amplitudes connected by 
(L + 1) cut legs. At six loops, we consider a seven-particle cut in the s = (k\ + k 2 ) 2 channel. 
The seven cut legs carry momenta p$,P6, . . . ,pu. The six dual loop momenta 2 5 , 2 6 , . . . , x w 
are then defined by, 

X45 = P5 , ^56 = P6 , X e7 = p 7 , 2 78 = p 8 , 2 89 = p g , Xg jl0 = p W . (2.6) 

The key dual conformal properties follow from the behavior of the integrand under dual 
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FIG. 2: Graphs 36 through 68 for the planar six-loop four-point amplitude. 



coordinate inversion, which maps 
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In four dimensions, dual conformal invariance requires that each term in the integrand scales 
as 
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The integrands of planar MSYM in D dimensions have been shown to transform in exactly 
the same fashion to all loop orders, at least for D < 6 37], |38[. This property is sufficient 
for our purposes, since we are mainly interested in the integrand in D = 5. 

The (L + l)-particle cuts can also be used to generate the complete list of graphs needed 
at six loops. One considers all possible sewings of two tree-level cubic graphs that appear in 
these cuts [l2|. (We modify the procedure slightly compared to ref. 12| by including only 
diagrams with cubic vertices.) In principle there are dual conformal graphs with four- or 
higher-point vertices that are not generated by the product of tree graphs of the (L + 1)- 
particle cuts; however, all such potential contributions, including those not detectable in the 
(L + l)-particle cuts, can be assigned to graphs with only cubic vertices by multiplying and 
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FIG. 3: The six-loop cut evaluated in D = 6 in ref. 37]. 



dividing by appropriate propagators. The construction of the potential numerators of each 
graph is then accomplished conveniently using the dual-graph representation, which exposes 
the dual conformal properties more simply. 

Given a dual graph in the sewing, we construct the possible monomials M^- as products 
of dual- momentum invariants xfj. We keep only those terms with the dual conformal 
scaling dictated by eq. (12. 8p . To determine the rational-number coefficients we use gener- 
alized unitarity. A large number of coefficients are easy to identify, essentially by inspection, 
because the corresponding unitarity cuts are so simple. In particular, contributions with 
either a two-particle cut or a box subdiagram can be written down immediately, following 
the discussion in refs. [!, 11 



In addition, many of the coefficients vanish. All but one of the vanishings can be 



identified using the observation of ref. |34j that when the external momenta are taken off-shell 
the integrals must be infrared-finite in four dimensions. The sole integral with a vanishing 
coefficient which cannot be identified in this way is the integral displayed in Fig. 12 of 



ref. [121 ] - It consists of two identical three-loop three-point integrals, containing only box 
subdiagrams, and connected to each other by one common external leg. 

The unitarity cuts include a sum over states in the supermultiplet for each cut leg. For 
generic cuts, the state sums are straightforward to implement numerically in four or six 



dimensions |37|, |39|, |40(. We have evaluated all four-dimensional cuts that decompose the 
amplitude into a sum of products of three-, four- or five-point amplitudes, as well as a 
variety of cuts involving six-point amplitudes. These cuts suffice to uniquely determine 



all the ciij. The two-particle cut |3[ and "box cut" |ll| are valid in D dimensions; hence 



all contributions to the integrand that are visible in such cuts are valid in any number of 
dimensions. In addition, a rather non-trivial cut of our expression, shown in fig. El was 
computed previously [13] using the six- dimensional spinor-helicity formalism of ref. (41| and 
superspace of ref. [42]. Thus this cut is valid for D < 6. This provides a highly nontrivial 
confirmation that the integrands are valid in D = 5. 

We have performed a variety of consistency checks on the integrand. The unitarity 
cuts offer highly nontrivial self-consistency checks, because the same monomial My can be 
visible in multiple cuts. As already mentioned we compared our integrand result against 
that of ref. [io| . As a further non-trivial check we confirmed that the unphysical singularities 
described in ref. [43[ all cancel. 

Some of the numerators iVj are quite simple. For example, for the graphs labeled (1) and 
(2) in fig. HI they are just 
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(2.9) 



Other numerators are more complex. For example, the numerator of graph (62) with the 
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FIG. 4: A few sample graphs with labels corresponding to the labels in eqs. (|2.9|) and (|2.10p and 
in the ancillary file [44J. The external momenta are outgoing. 



labels in fig. H] is 

^62 = ^Stplpl + t{k 1 +p 21 )\k 2 -p 18 )Y 9 pl 7 -tpl 7 p 2 20 {k 1 +p 21 ) 2 {k 3 + Pl3 ) 2 (2.10) 

- splpKh +Pn) 2 (h +P23) 2 ~ spjjploih + p l7 ) 2 (k 3 + p 23 ) 2 

- stp^h + p 2 i) 2 (p 20 - pg) 2 + t(ki + P2i) i (h + pi 3 ) 2 {p 2 o - p 9 ) 2 
+ s(fci + pnfih + p 21 ) 2 (k 3 + p 23 ) 2 (p 20 - P9) 2 ■ 

In these expressions we have chosen labels that line up with the ones in the ancillary text 
file [HJ]. The symmetry factors of these graphs are Si = 4, S 2 = 2 and Sq 2 = 1. The 
complete sets of diagrams, numerators Ni, and symmetry factors Si, are included in the 
ancillary file 44| . 



III. FROM THE AMPLITUDE TO VACUUM DIAGRAMS 



Because the UV divergences arise from integration regions in which the loop momenta 
are parametrically much larger than the external momenta, extracting the UV divergences is 
a much simpler task than integrating the complete amplitude. We follow the same strategy 



as m our previous papers 



10, 23 



based on Taylor expanding the integrands in small 
external momenta and then integrating the resulting vacuum integrals [21], 22| . The present 
case is relatively straightforward to analyze in the sense that the six-loop amplitude contains 
no sub divergences in D = 5, and because the expected overall divergence is manifestly log- 
arithmic. However, the high loop order makes the loop integration for the vacuum integrals 
highly nontrivial. 



A. Obtaining the vacuum diagrams 

We consider the individual integrands Xj appearing in eq. f)2.3p in the limit of small 
external momenta kj, j = 1, 2, 3, 4. We let kj — > ekj, and then expand in the small parameter 
e, keeping only the leading order. This reduces each Xj to a sum over six distinct vacuum 
integrands, 

l^ek^pt) e 2 (sA l!X + tB hX )V {x \ Pl ) + 0(e 3 ), (3.1) 

xg{a,b,c,d,e,f} 
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FIG. 5: The six distinct vacuum diagrams that appear in eqs. (|3.3p and (|3.4p . Each dot indicates 
that the corresponding propagator should be squared (doubled) in the integrand. The five "tensor" 
integrals have numerator factors that are indicated by the prefactors. The numerator factors are 
built from momentum invariants involving a subset of the loop momenta, labeled by h,l2,h,h- 



where Ai )X and B^ x are rational numbers determined by the expansion. (We will not list 
these coefficients separately for each diagram.) After the above vacuum integrands are 
integrated over the six loop momenta p^,Pq, . . .pio in D = 5 — 2e, with the measure 
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5-2e 



Pi 



(2tt)' 
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we obtain six vacuum integrals, V^°\ . . . , shown in fig. [5j These integrals have 
numerator factors which are indicated to the left of each graph, and either one or two 
doubled propagator, whose location is indicated by a dot. The integrals V^ x > contain no 
subdivergences; each integral has a single overall UV divergence in D = 5 when all six loop 
momenta become large. Hence the integrals have only simple poles in e. 

Collecting the contributions from the 68 distinct integrals in the six-loop amplitude 
eq. ( 12. 2p . we obtain the following UV divergence 



.4 



(6) 



D=5, div. 



= 6stuAf*(l, 2, 3, 4)(^ (a) + V {h) + 2V {c) + 4V {d) + 2V {c) - 2V {1) ) . (3.3) 



10 



FIG. 6: The canonical vacuum integral. All contributions can be expressed in terms of this diagram. 
The differences between each contribution can be assigned to the one-loop subintegral indicated 
by the shaded region. 



This simple expression for the divergence appears to be nonzero at first glance. However, 
we should expect that these six integrals are not independent so there may be cancellations 
among them. While all six integrals are positive definite after Wick rotation, enters 
with a coefficient of the opposite sign from the others. Thus no conclusion can be reached 
as to whether this expression vanishes or not without a careful analysis. 

We begin the analysis by collecting all integrals under one common integration. By 
multiplying and dividing by propagators, all contributions can be brought to the form of the 
single vacuum graph displayed in fig. |6j Then by appropriately relabeling each term from 
the original Pi to another set of loop momenta, all contributions can be rearranged to differ 
in only a single one-loop subintegral, indicated by the shaded region in fig. [61 The momenta 
h,l2,h,U are carried by the indicated one-loop subintegral. The momenta m 1 ,m 2 ,m 3 ,m4 
are the four momenta external to the one-loop box subintegral. In addition, we need two 
other independent loop momenta, which we take to be q\ and q2, as indicated in fig. [6j Of 
these momenta, we take l±, m 2 , m 3 , m 4 , q±, q 2 to be the six independent ones. In order to 
make the analysis easier, we may symmetrize the contributions over the automorphisms of 
the vacuum integral. 

After appropriate relabelings, we can express the sum over all vacuum integrals as a single 
integral, 



V = V {a) + V (b) + 2V {c) + 4V id) + 2V {c) - 2V r(f) 
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+ <?2) 2 (<?i - "^4) 2 (gi - m 4 - m 1 ) 2 (q 2 - m 2 - m 3 ) 2 (q 2 - m 3 ) 2 
where e is the dimensional-regularization parameter and the vacuum "numerator" is 
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A^vac = (Zi + h)' 
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using the momentum labels in fig. [6j that is, h = h — m 2 , h = m 2 + m 3 — h, U — —mi — l\ 
and WI4 = —mi — m>2 — m^. By a slight abuse of convention, we call this a numerator even 
though it is nonlocal. The numerator depends only on the momenta internal and external 
to the one- loop box sub diagram indicated by the shaded region in fig. El 

Because of the minus sign in the last term — corresponding to —2V^' in eq. (13. 3p - 
it is easy to see that this integrand is not positive definite, even after Wick rotation. For 
example, for l\ = —l 3 all terms but the last one vanish, making the integrand negative; 
and in the region where < m\ is much smaller than all other momentum invariants, the 
terms with \jm\ factors will dominate, making the integrand positive. Thus, even after 
combining all contributions into a single integrand, there does not appear to be a simple 
way to determine the positivity, or vanishing, of the integral. 



B. Simplifying the vacuum integrals via consistency relations 

To simplify the expression further we need to identify relations between the different 
vacuum integrals. The integral identities that we need are related to integration-by-parts 
identities (IBP) [24] ]. although they are only valid for the leading 1/e UV pole in the critical 



dimension, D = 5 — 2e. These consistency relations are obtained by demanding that 
different loop-momentum parameterizations of the integrals lead to the same final results. 
Using these relations we can both simplify the UV divergence, and give important cross- 
checks of the numerical evaluation. The latter use is particularly important, because the 
relations give independent estimates of the numerical uncertainties. The integral consistency 
relations also offer a potential path to finding a positive-definite expression for the divergence. 

We now sketch a derivation of the consistency relations. For each of the 68 graphs 
describing the amplitude, we have simple relations from the shift invariance of the integrals, 

for each q m £ q 6 }, and Xj, for i — 1, . . . , 68, are the integrands of the distinct graphs. 

The Ni are general numerator polynomials in the momentum invariants of the graphs. These 
polynomials are chosen to generate useful identities and are not the numerators iVj of the 
amplitude. A judicious choice of a large set of Ni will lead to a large set of linear relations 
for various vacuum integrals, which will include (but will not be limited to) the desired 
integrals V^ x \ Below we describe such a judicious choice of numerators. 

The identity (13. 6p follows because the q\ momentum dependence of the integrands is 
completely removed by a change of variables in the measure of eq. (I3.6p . In fact, eq. (13.61) is 
simply a statement that the integrals are reparametrization invariant under constant shifts. 

Next we expand the integrands in small external momenta, kj — > ekj, with e a small 
parameter. In doing so, it is convenient to treat the integrands as belonging to equivalence 
classes controlled by the reparametrization freedom, 

l i [N i ](ek j ,p l + q t ) ~ UNijiekj^i) . (3.7) 

Expanding the two sides of eq. (13. 7p in e would not yield any nontrivial equations, only trivial 
reparametrization relations for vacuum integrals. However by combining the reparametriza- 
tion freedom and the small momentum expansion, that is, by letting q\ = e^. cyfcj, we get 
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nontrivial relations between different vacuum integrals, 

li[Ni] (ekj,pi + Cl ^J ~ ^iWi]{ekj,pi) , (3.8) 

3 

where cy is an arbitrary integer-valued 6x3 matrix, where different choices will generate 
different relations between vacuum integrals. The expansions in e of the two different inte- 
grands will differ considerably. But by the reparametrization freedom, the two sides must be 
equivalent as integrands, or equal after integration. Therefore the various vacuum integrals 
that arise from integrating the coefficients of each element of cy in eq. (I3.8P must satisfy 
nontrivial consistency relations. 

It is important to make judicious choices for the integral numerators JVj used to generate 
useful consistency relations. For example, using the original JVj appearing in the integrand 
of the amplitude is not a good choice, because their divergences are manifestly logarithmic 
in D = 5 — 2e. The UV divergence of a logarithmically-divergent integral is given by the 
leading term in e. This term is always insensitive to the shift in the loop momenta qi, 
and so the available consistency relations become trivial. The first nontrivial relations are 
obtained using numerators iVj containing one additional power of a loop momentum, which 
give linearly-divergent integrals in D = 5 — 2e. The relevant vacuum integrals relations are 
obtained from the next-to- leading term in e, which will differ on both sides of eq. (13 .7p . 
Numerators iVj that give rise to quadratic divergences are also useful for extracting integral 
relations. However, numerators with an even higher degree are less helpful, since after taking 
derivatives with respect to e they give rise to vacuum diagrams with three or more doubled 
propagators. These are outside the class of integrals that we are interested in; the vacuum 
diagrams in fig. [5] have at most two doubled propagators. 

Furthermore, one should not choose Ni that give rise to sub divergences in D = 5 — 2e, 
because then the consistency relations may be contaminated by relations that are only valid 
for overlapping leading 1 je n UV poles. Specifically, for integration in D = 5 — 2e dimensions, 
UV sub divergences are possible in principle for two- and four-loop subdiagrams. Any JVj 
generating such a subdivergence should be eliminated from the set of choices, because it will 
not produce any useful identities. 




FIG. 7: Two simpler vacuum integrals that appear in the UV divergence after using integral 
identities. The third integral that enters the divergence is fig. EJf). 

Generating a sufficient set of consistency relations then comes down to varying the iVj 
polynomials for an appropriately large function space, without exceeding available compu- 
tational resources. This includes varying the matrix Qj that controls the reparametrization 
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of the integrand, and then, as explained, expanding the integrals in small external momenta 
and demanding that the expansion is consistent for different choices of Qj. 

After generating about 1000 independent consistency relations, we found a much simpler 
three-term expression for the UV divergence of the planar six-loop four-point amplitudes. 
We then generated an additional 7000 consistency relations, and no further improvement was 
found. Thus a search for identities beyond the ones we found would probably be unfruitful, 
though we have not proven that they do not exist. 






(k) 




FIG. 8: Six scalar vacuum integrals that appear at intermediate steps in integral identities that 
yield a simplified UV divergence. 



Using the derived consistency relations, the five tensor integrals . . . , V^ e ' can be 
reduced into eight scalar integrals with no loop momentum in the numerator, plus one 
relatively simple integral with only two powers of loop momentum in the numerator. One 
of the scalar integrals is V^' from fig. The eight new integrals are shown in fig. [7] and 
fig. [BJ The most useful of the derived consistency relations are: 



y(e 



2U (f) + 2U (g) - 4U (i) + 2U (k) - 2U (1) , 
y(f) + 3y(s) _ 4V (i) + 2U (k) - 2U (1) , 

Iy(f) _ Iy(h) + yd) _ 2y (m) + y(n) 

2 2 

ly-(f) + 2V® - U (k) + U (1) 
2 

_yG) + 2U (m) - V (n) . 



(3.9) 



After applying the integral identities (13. 9p to eq. (I3.3p . we obtain the following simplified 
form for the UV divergence, 



.4 



(6) 



D=5, div. 



QstuAT c (l, 2, 3, 4) (10 V {1) + 5 V {s) - V {h) ) 



(3.10) 



Unfortunately, the coefficient of vacuum integral has a relative negative sign, so this 
simplified form is also not positive definite. One may wonder if there exists a different choice 
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FIG. 9: An extra scalar diagram used for checking integral identities. 

of basis of vacuum integrals that that allows the divergence to be expressed in a positive- 
definite form. While this possibility cannot be excluded by our analysis, we have not found 
such a representation, and it is quite likely that eq. ( 13.1 Op is the simplest integral basis 
decomposition. 

In addition to the needed relations, and as a cross check of the numerical evaluation in 
section |V] we offer one very simple integral identity between three scalar integrals 

V {o) = ^V® + V {i \ (3.11) 

where is shown in fig. EJ 



IV. REVIEW OF SECTOR DECOMPOSITION 

Given that even the simplified form (13.101) still leaves open the question of whether the 
amplitude diverges in D = 5, and because analytic techniques are not yet powerful enough 
to cope with generic six-loop integrals, we have resorted to a numerical evaluation of the 
relevant vacuum integral using the method of sector decomposition. In this section we review 
this method, focusing on the salient features needed in our calculation. 

If the amplitude under consideration is UV divergent, numerical evaluation with even a 
modest (say 5%) accuracy will suffice to prove it beyond any doubt. Indeed, here we shall 
provide such a numerical proof for MSYM in D = 5, settling the question of the potential 
finiteness of this theory. Even had it turned out finite, numerical analysis can provide 
important evidence in favor of this hypothesis. We are interested in the leading-logarithmic 
divergence of an integral with no sub divergences. After Feynman parametrization of the 15 
propagators and integrating over the loop momentum, the overall divergence appears as a 
coefficient of a convergent 14-dimensional parametric integral. Eqs. (13. 3 p or ( I3.10p are fairly 
simple and it might seem to be an easy job to estimate it by Monte Carlo methods. If we 
estimate the integral as an average over N uniformly distributed samples, we might imagine 
that the statistical uncertainty of such an estimate would be a/\/N, where a is the standard 
deviation of the integrand (a 2 = J d D x(f(x) — f) 2 with / the average value). An accuracy 
of 1% would seem easily attainable. 

However, the situation is more complex. The problem is that Feynman integrals in 
general, including eqs. (13 ,4p and (I3.10p . are not sufficiently convergent because of endpoint 
singularities. UV divergences in themselves are not a major problem, and are usually dealt 
with by subtracting a simpler integrand with the same divergent behavior. In the case at 
hand, even this is not necessary because the coefficient of the divergence is given by an 
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(absolutely) convergent integral. However, the square of this integrand is not integrable. 
Because of this, the variance a 1 will diverge, and Monte Carlo estimation with the uniform 
measure has uncontrolled errors. In practice, it does not work. 



A. Sector decomposition overview 



A straightforward way to deal with this problem is to carry out Monte Carlo integration 
with a sampling measure designed to overcome this problem. Let dfi(x) be the sampling 
measure. Then we rewrite 



d D x f(x) 



d{i(x) 



d D fj,(x) 



det 



d D 



x 



/(*) 



= / dn{x) F(x) 



(4.1) 



where F(x) combines the original function f(x) with a change of measure factor. A particular 
case would be to take a new set of variables y functionally related to x, and d\i to be uniform 
measure d D y, in which case this factor would be the Jacobian for the change of variables. 
The error estimate for a Monte Carlo integral with this sampling measure is ap/y/N where 
op is the new standard deviation, 



J dii{x) {F(x)-F) 2 , 



(4.2) 



where F is the average value of F. By reducing the variance, one can improve the accuracy 
of the Monte Carlo. In principle one could change variables to absorb all of the variation 
of f(x) into the measure to make F(x) constant, in which case the exact result would be 
obtained from a single sample. Of course, in practice such a change of variables would be 
prohibitively difficult. 

For an integrand with power-like singularities, such as 



dxx a ~ 1 (l + f(x)) 



(4.3) 



we could use the change of variables y = x a /a, which absorbs the singularity into the 
measure and does not much complicate the integrand. The problem with applying this to a 
multi-dimensional Feynman integrand is that it has many different power-like singularities, 
arising from the many orders in which one can take the different parameters Xi to zero. 

The solution is to decompose the integration region into subregions or "sectors," each of 
which has at most one singular behavior of this type. If we can do this, we can apply the 
change of variable eq. (14. 3p . or its multivariate generalization, to regularize the integral in 
each sector. While this approach was long used in formal proofs of perturbative renormal- 
izability [45 44 7]. it seems to have been first used in numerical integration by Binoth and 
Heinrich [25J. Since then has been implemented in several computer pa ckag es for numeri- 
cally evaluating Feynman diagrams, starting with Bogner and Weinzierl [26| and including 
refs. (27l-[29|. Let us explain the basic ideas, leaving the details specific to our computation 
to section IVl 
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We begin with an elementary example (from section 2 of ref. [48|), the two-dimensional 
integral 

The form of the denominator makes the limit x, y — > hard to control. Although in this 
simple example we could change variables to (say) x and w = x + (1 — x)y, this option will 
not be available for more complicated integrands. 

Rather, we split the integration region into two parts, region 1 with x > y and region 2 
with y > x. In region 1, we can make the change of variables 

x — x' , y = x't' , (4.5) 

turning the integration region into < x',t' < 1. Similarly, in region 2 we take 

x = y't' , y = y', (4.6) 

again turning the integration region into a square. The integral becomes (suppressing the 
primes) 

1 rl x Q+ ^t /3 f l f 1 t a V a+l3 

dx / dt- + dy dt- f- -. (4.7) 



/o Jo l + (l~x)t J J -y)t 

Now the nature of the singularity is manifest in the leading monomial terms, because the 
complicated denominator goes to 1 in the singular region. 

The same idea can be applied to a function of N variables, call these Xi with i e [1, N]. 
The integration region is decomposed into N subregions labeled by a G [l,N] and defined 
by the inequality 

x a > Xi, i ^ a. (4.8) 

In the a'th sector we redefine 

Xi = Xax'i , a ^ i , (4.9) 

to turn the subregion back into a unit cube. This will allow pulling out an overall singular 
behavior controlled by x a . Of course, the resulting integrand might still have a complicated 
singularity in the other variables. This must be dealt with by iterating the procedure and 
dividing the subregion into further subregions. Mathematically, this operation is called 
"blowing up" the singularity. 

In these subsequent subdivisions, one need not include all N variables in the blowup; 
one can instead take a subset of the variables and apply the same procedure. These choices 
might be used to simplify the result, or might even be needed in order for the procedure to 
terminate. The hope is that by choosing an appropriate sequence of these operations, one 
can find a finite set of subregions in each of which the integrand takes a simple form, as in 
eq. (H2D. 

This procedure may be familiar to some readers from its use in algebraic geometry and 
string compactification, and the following remarks are addressed to them. In complex alge- 
braic geometry, one can blow up an arbitrary point p in an iV-dimensional space, replacing 
it with a CP^" 1 . This is done by taking coordinates in which Xi(p) = and applying the 
same changes of variables; the a'th subregion corresponds to the coordinate patch on CP^ 1 
in which we can take x a = 0. 

Suppose that the integrand is a rational function with denominator D(x). The singularity 
is then the set of all points satisfying D(x) = 0. Let us denote this singular set as A. Since 
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for a Feynman integrand the function D{x) is a polynomial, the set A is by definition an 
algebraic variety, meaning the set of solutions of a system of polynomial equations. In fact it 
is a hypersurface, since we are setting a single polynomial to zero. In this context, a natural 
thing to try is to blow up the space containing A, to a variety n : X — > C N so that 
7r _1 (A) is nonsingular, meaning that it is defined locally by a single equation / = with 
df 7^ 0. If we can do this, then the singular behavior of the integrand will simply be 1/f 
(or perhaps some power of this). By taking / to be a local coordinate, we would accomplish 
our goal of realizing the singular behavior in a particularly simple form. 

For present purposes, the main result of this mathematics is Hironaka's theorem on 
resolution of singularities, which states that any singularity of an algebraic variety can be 
resolved (made nonsingular) by a succession of blowups. Furthermore, there are algorithms 
for concretely finding the resolution. Thus, we can use a blowup algorithm to resolve the 
singular locus A, providing; the multiparameter generalization of eq. 04.7p . 

While in ref. 26| this idea was used to give blowup algorithms which are guaranteed to 
terminate, these algorithms tend to produce a number of subregions which are exponential in 
the number of variables — this is perhaps natural as the number of orderings of N variables 
is N\. This number of subregions would be computationally infeasible for N = 15. Bogner 
and Weinzierl [26| also proposed a simpler heuristic algorithm which, while not guaranteed 
to terminate, produces a simpler solution when it does. Another heuristic algorithm was 
proposed by Smirnov and Tentyukov [28|, which we now describe. 



B. Heuristic sector decomposition 

We recall from textbooks (e.g. ref. [49j] ) that the denominator of the Feynman integrand 
for a vacuum integral is the D/2 power of the Kirchoff polynomial of the graph, 

t ifr 

where the sum is taken over the spanning trees T of the graph T. We are interested in the 
limiting behavior of U as combinations of the variables go to zero. This behavior is encoded 
in its Newton polytope. Let deg be the degree of a monomial, considered as a vector in Z N , 
so that 

dega:™ 1 ^ 2 ...x n N N = (n 1; n 2 , . . . ,n N ) . (4.11) 

The Newton polytope of U is the convex hull of the degrees of each of its terms; in other 
words, it is the set of all points in R which can be obtained as linear combinations of these 
degrees with non-negative coefficients. 

We need to desingularize each limit which takes a subset of the variables to zero. We 
now assume that U is a polynomial with no constant term, so that every monomial in U 
will go to zero for some such limits. However, many of the monomials are subleading and 
do not control any limit: if a monomial Mi is the product of another monomial M 2 with a 
monomial of non-negative degree, it is subleading. In terms of the degrees, this requires 

deg M x - degM 2 > 0, (4.12) 

for every component. 
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We refer to the points deg Mi which do not satisfy eq. (I4.12p for any deg M 2 7^ deg Mi as 
the "low points" of the polytope. If the polytope has a single low point, then by factoring 
out the corresponding monomial, one obtains a polynomial with a nonzero constant term. 

If there are multiple low points, we need to do a blowup. A blowup on a subset S of the 
variables subdivides the current sector into \S\ sectors. In the a'th sector one applies the 
change of variables (14 .9|) . This change of variables operates on the Newton polytope as 

v v + e a (xs - e a , v) , (4.13) 

where xs is the vector whose components are 1 for i G S and for i ^ S, and (v,w) = 
Y2i v i w i- We will then be able to factor out a common monomial; in other words shift the 
entire polytope in a way that keeps it in the upper quadrant. If the result includes the 
origin, we are done with this sector; otherwise we apply the same procedure recursively. 

The next problem is to decide which subset of variables to involve in the blowup. The 
goal is to eliminate as many low points as possible; however it is better to leave the variables 
which do not contribute to this goal out of the blowup, in order to maximize the degree of 
the monomial that can be factored out. A condition on the subset S which favors this is to 
project the polytope onto the subspace V5 and require that the low points of the projected 
polytope linearly span this subspace. The heuristic algorithm is simply to choose, from the 
subsets satisfying this condition, the subset S that has the largest value of Yli^s 

This heuristic algorithm is too simple to desingularize general polynomials, including 



examples that are given in ref. |26j. If applied to these examples, it will go into an infinite 
loop. 1 However, the heuristic algorithm works for the class of polynomials in eq. (I4.10p ; that 
is, scalar vacuum integrals with no doubled propagators. It also turns out to work for vacuum 
diagrams with no IR divergences and simple numerator factors, including the integrand (13.41) . 
As observed in ref. jHo|, it works because it reproduces the results of a canonical sector 



decomposition procedure, which associates sectors with maximal forests [46|, |47|, l5jj . While 
we leave the details for the references, a maximal forest is a hierarchical decomposition of 
the graph into a set of subgraphs satisfying certain conditions (each pair of subgraphs (7, 7') 
must obey one of the relations 7 C 7', 7' C 7 or 7 fl 7' = and 7 U 7' can be disconnected 
by removing a single vertex). It can be shown that a maximal forest for a diagram with L 
loops and E edges contains E trivial subgraphs (single lines) and L nontrivial subgraphs, 
and the associated sector involves L blowups each on distinct variables. 

If the heuristic algorithm is reproducing this decomposition, then since every sector in- 
volves a succession of L blowups on distinct variables, the algorithm is guaranteed to ter- 
minate with at most N\/(N — L)\ sectors. For N = 15 and L = 6 this is 3603600 which is 
not much larger than the actual numbers we obtained. For N = 18 and L = 7 it is about 
1.6 x 10 8 . 



1 The "Strategy X" of ref. [26( and "Strategy S" of ref. [29( are somewhat more sophisticated and can handle 
these cases. 
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V. NUMERICAL RESULTS AND COMPUTATIONAL DETAILS 



A. Numerical integration 

There are several software packages for carrying out sector decomposition and subtraction 



of divergences, and for i nteg rating the resulting expressions numericall y |27H 29|. We used 



the package FIESTA 2 [29|, written in a combination of Mathematica [52[ and C++ [53 
and which can take advantage of multiple processors. However, a six-loop diagram is too 
complicated to evaluate directly. For example, the sector decomposition (which is not very 
parallelized) takes about a day per primary sector to run in Mathematica on a modern 
computer, and produces a total code for the integrands which takes up hundreds of gigabytes. 
Thus the computation must be split into smaller parts to make it feasible. 

To do so, we only need a small portion of the FIESTA 2 software, and we have extracted 
this portion and adapted it to our purpose by hand. We did the sector decomposition on 
a small (10 node) cluster, and then performed the numerical integrations on a large (1000 
node) cluster, using the adaptive quasi-Monte Carlo integrator Vegas [Hj]. An important 
element in FIESTA 2 is the CIntegrate package (see ref. [2^|, appendix F), which accepts 
a symbolic algebraic expression of the sort that can be produced easily by Mathematica, 
and compiles it into pseudocode that can be executed efficiently in C. Using this package, 
we were able to break down the computation by having Mathematica produce a file for each 
sector containing its integrand, which could be passed to a VEGAS integration program. 
The integrals are of course completely independent, so this step is easy to parallelize. This 
approach also allowed us to keep the partial results for every sector, which helped in de- 
bugging and uncertainty analysis. The main cost was the need for 1-2 TB of disk storage, 
which is not large these days. 

A rough estimate of the total running time can be obtained by multiplying the number of 
sectors (about 10 6 here) by the time to integrate a sector, divided by the number of nodes. 
With a regular (low variance) integrand, the uncertainty as a function of the number N of 
samples goes somewhere between N~ l l 2 for Monte Carlo and iV -1 for quasi-Monte Carlo in 
low dimensions. Vegas did not need more than 50, 000 samples to achieve our requested 
relative precision of 10 -3 in any sector, and this took 2-3 seconds to do. Thus a million 
sector integration took less than an hour on the cluster. 

There were a number of reasons that this success was not guaranteed from the start. 
Even once we knew that we had of order 10 6 sectors, the next possible pitfall was that the 
integral might be small due to cancellations between larger results in individual sectors. 
The relative uncertainty of course depends strongly on the relative signs of the intermediate 
results; in the best case (a single sign) we might hope to gain a further statistical 1 / ^/N scctors , 
while in the worst case the result might be comparable to the largest statistical uncertainty 
in a single sector (which could of course be improved by taking more samples) or even an 
uncontrolled systematic uncertainty. In fact, it appears that, at least in our computations, 
sector decomposition does lead to significant systematic uncertainties, as we will see below. 

The actual situation is best judged by looking at the partial results, which we graph for 
diagram U (g) in fig. [TOJ There are 1224600 sectors, all of which gave positive contributions. 
The same is manifestly true for and (since these integrands were positive), so there 
is no problem with cancellations. 

Although a few sectors give anomalously large values, the largest sector value is only 
about 40 times the average sector value, while the largest uncertainty is about 180 times the 



20 



Sector Value divided by Mean Sector Value 




FIG. 10: The sector integrals for diagram , in order of increasing numerical value. The dark 
bottom (black) line gives the sector value, normalized by dividing by the average sector value. The 
top light (red) line gives the natural logarithm of the normalized sector value. 
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FIG. 11: The total contribution of sectors 2™ through 2 n+1 , ranked in order of decreasing numerical 
value. 



overall standard deviation. One also sees that the distribution of sector integrals is roughly 
exponential except at the ends. This point is made in a different way in fig. (TTJ In this 
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plot the sectors are listed in decreasing numerical order, with the graph showing the total 
contribution of each subset which has an exponentially growing number of sectors. Because 
the totals in each subset are roughly constant, this shows that the sector contributions 
decrease exponentially. Thus, while one cannot reduce the calculation to evaluating a smaller 
number of dominant sectors, there is structure which might be exploited to speed up the 
computation. For example, if one simply fits the distribution to an exponential and integrates 
that, one gets within 20% of the actual total. 

All this seems (at least naively) consistent with the uncertainties being entirely statistical, 
in which case it would be appropriate to add them in quadrature to get the statistical 
uncertainty estimate for integral of 6 x 10~ 4 %. However, in previous use of sector 
decomposition, it has been observed that a purely statistical combination underestimates 
the true uncertainty. A more conservative estimate would be to treat the uncertainties as 
100% correlated — what we will refer to as a "systematic" uncertainty estimate. In this 
case, we should add the uncertainties in every sector. For integral V^ 9 \ this error is 0.29%. 

In Table HI we collect our results for the various vacuum integrals defined in section IHH 
along with uncertainties that have been estimated by assuming 100% correlation among 
sectors. These results can be used to confirm the consistency relations (j3.9P ; the degree of 
agreement is in turn a cross check of the integration uncertainties. 

For example, consider the first equation in eq. (13. 9p . This relation evaluates to the 
numerical values 

1.3958 ~ 2 x 0.7631 + 2 x 0.8183 - 4 x 0.5967 + 2 x 1.1493 - 2 x 0.8391 = 1.3964, (5.1) 

which is certainly acceptable and consistent with the systematic uncertainty hypothesis. The 
other consistency relations in eq. (13. 9 p work to a similar accuracy. Thus the combination of 
sector decomposition with the Vegas adaptive numerical integration appears to introduce 
systematic uncertainty, in the sense that the error is correlated between different sectors. 
Presumably this has to do with the modeling of the integrand at the sector boundaries. 

Taking the uncertainties as 100% correlated and conservatively adding them directly 
(instead of in quadrature), the final result for the UV divergence is 



(6) 



D=5, div. 



X6(lOV {i) + 5V {S) -V^ 

10 x (0.7631 ± 0.0015) + 5 x (0.8183 ± 0.0024) 
(0.2762 ± 0.0008) 



1 X 6 



£(4*; 
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where 



7(4^ (68 ' 68±0 - 17) ' (5 - 2) 



X = -stuAl ree (l, 2,3,4) . (5.3) 



It is clear from this result that the coefficient of the UV divergence is nonzero, well within 
the integration uncertainty. This proves that MSYM is perturbatively ultraviolet divergent 
in D = 5. 
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Integral Value Uncertainty (Value) / 



y w 


1.3958 


0.0043 


5.05 


T r(W\ 


1.4522 


0.0079 


5.26 


rr(c\ 

yv~i 


1.3346 


0.0069 


4.83 


y\ u ) 


1.2643 


0.0029 


4.58 


trie) 


1.1935 


0.0026 


4.32 


yv-) 


0.7631 


0.0015 


2.76 




0.8183 


0.0024 


2.96 


rr(h) 

y\ u ) 


0.2762 


0.0008 


1 


yv) 


0.5967 


0.0012 


2.16 


yG) 


1.1490 


0.0020 


4.16 


^(k) 


1.1493 


0.0019 


4.16 




0.8391 


0.0015 


3.04 


y( m ) 


1.8600 


0.0028 


6.74 




1.3755 


0.0022 


4.98 




0.9790 


0.0017 


3.55 



TABLE I: The data used to numerically verify the integral consistency relations, as a means of 
assessing uncertainties in the numerical integration. The columns labeled by "Value" and "Uncer- 
tainty" are multiplied by e(47r) 15 . 



B. Extrapolating in loop order 

Maximally supersymmetric Yang-Mills four-point amplitudes have a smooth analytic be- 
havior as a function of dimension, through at least six loops. As seen from their explicit 
forms the only dependence on the space-time dimension is in the loop momentum integra- 
tion measure. Might this property lead somehow to a simple functional form describing the 
numerical values of the divergences in the critical dimensions, as a function of the number 
of loops L? To check this hypothesis, we plot the known values of the divergences in fig. [12] 
represented by the dots. In doing so we extract some simple overall factors, defining the 
numerical constant pi by 

A™\ = H ~ l)L ~ 1X p, (5 4) 

^4 lz) = 4+6/L,div. e (4^)2X4-3 ^ L ' 

where X is defined in eq. ( 15. 3p . 

The L = 1 value is not plotted because it does not obey the critical dimension bound 
given in eq. (11.11) . and its divergence in D = 8 differs kinematically from eq. (15 .4p by a factor 
of 1/u: 

\ X \ \ 

4 'a=8.div. = e(4^6«- (5 ' 5) 
The values of the divergences in.D = 4 + 6/L from two through five loops have been given 
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FIG. 12: The dots indicated the numerical coefficients in eq. (|5.6p of the UV divergences in the 
critical dimensions. The solid (blue) line is the result of fitting the parametric form in eq. (|5.7[) to 
the displayed results for L = 2, 3, 4, 5, 6. The dashed (purple) line is a fit to the parametric form in 
eq. (|5.1U|) . The lower panel shows that the relative error between the points and the fit in eq. (|5.7|) 
is within 1%. 



previously (H 0, [ljj, 13]. Correcting a couple of overall signs, they are 



P2 


71 

20' 




1 




P3 


3' 




= 6 


"512 


Pi 


5 


P5 


~ 9.537, 



'3\4 



^r(i) 3 r ( |)r ( i; 



1.553 



(5.6) 



where the expressions through L = 4 are exact. The L = 5 expression is approximate, but 
it is accurate to the digits given. 

The linear behavior beyond L = 2 in the upper panel of fig. [T5] make clear that the 
coefficients of the divergences have an approximately exponential behavior. This observation 
motivates a simple Ansatz for the approximate form of the divergences at any loop order 
L > 2, 

PL ^6 lC f +ai/L . (5.7) 
The solid curve in the upper panel of fig. [12] is based on eq. (15.71) with the parameters 



^ = 3.99 



&! = 1.74 x 10" 



ci = 9.77. 



(5.8) 
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Interestingly, a nearly equally good fit is given by the following analytic form, which contains 
remarkably simple constants, 

PL * (n 2 ) L+A/L ~ en/2 • (5-9) 

Since we do not know the precise functional form, for the purposes of extrapolating to 
higher-loop orders, it is useful to compare this to a different Ansatz, 

p L ^a 2 + b 2 c%, (5.10) 

where again a 2 , b 2 and c 2 are parameters. The dashed curve in fig. [12] corresponds to 
eq. ( 15.1 Op with the parameters 

a 2 = 0.179, b 2 = 4.52 x 10" 4 , c 2 = 7.30 . (5.11) 

Before extrapolating to higher loops, an interesting exercise is to use eqs. (15. 7\i and (I5.10p 
to see how well they predict (or rather postdict) the obtained UV divergences for L = 5,6. 
Because the Ansatze involves three parameters, we need to use the values for L = 2, 3, 4 to 
fix the function. For the anstaz in eq. (I5.7P we obtain 

ai=4.06, h = 1.34 x 10" 5 , a = 10.2 . (5.12) 
Plugging this solution into eq. (I5.10p gives good predictions for the L = 5, 6 cases, 

p 5 ~9.91, p 6 ~74.0. (5.13) 
Similarly, for the Ansatz in eq. (I5.10p . we obtain 

a 2 = 0.127, b 2 = 6.22 x 10" 4 , c 2 = 6.92 , (5.14) 

and 

ps-9.99, p 6 ~68.4. (5.15) 

Presumably, the surprisingly good agreement between these approximate values and the 
calculated ones in eqs. (15.21) and (I5.6P is somewhat accidental. Nevertheless it does illustrate 
the remarkably good predictive power of this simple extrapolation. 

Using the fit parameters based on explicit results through six loops, we can easily predict 
approximate values for higher loops. For example, through L = 9 from the Ansatz (15. 7p we 
have, 

p 7 ~ 540 , p 8 ~ 4490 , p 9 ~ 38700 , (5.16) 

while the Ansatz (15.101) gives 

p 7 ~ 500 , p 8 ~ 3650 , p 9 ~ 26700 . (5.17) 

The small numerical integration uncertainty from our L = 6 result feeds into this fit, prop- 
agating a few percent spread in the estimates. Of course, the functional forms may be too 
naive, but the two different fits gives an indication of the spread in predictions for such 
extrapolations. 

Another interesting numerical observation from refs. 10, 13[ is that for an SU(N C ) gauge 
group, the ratio of the 1/iV^-suppressed subleading-color contributions to the leading color 
ones are fairly constant for L = 3,4,5 and takes a value of about 45. This observation 
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immediately gives us a prediction for the value of the divergence for the fully color-dressed 
amplitude (including nonplanar contributions), 



A {6) 



- ~ Tir^ 9 1A stAT e iV c 4 (68.68 N 2 C + 3100) 
div e (4tt) 15 4cv c i 



(5.18) 



X 



s(Tri 324 + Tr 1423 ) + t(Tr 1243 + Tri 342 ) + u(Tr 1234 + Tri 432 ) 



where Tr 4234 = Tr [T ai T a2 T a3 T a4 ] and we assume that only the leading-color and 1/N%- 
suppressed single-trace terms contribute, as is the case for L = 3,4,5. 

It is interesting to note that, at least through four loops, the divergences of M = 8 
supergravity are in the same critical dimension D c = 4+6/ L, and they are proportional to the 



same linear combination of vacuum integrals as the 1/A^-suppressed terms of MSYM [10l.l55 



If this link between gravity and gauge theories were to persist to all loop orders, then the 
critical-dimension agreement alone would imply the four- dimensional ultraviolet finiteness 
of the theory. It would be very interesting to directly check these ultraviolet divergence 
patterns in both gauge and gravity theories at as high a loop order as possible, in order to 
see if they could give insight into the UV properties of Af = 8 supergravity, and into the 
precise values of the ultraviolet divergences in MSYM in D = 4 + 6/L, to all loop orders. 



VI. CONCLUSIONS AND OUTLOOK 

We have shown that planar MSYM diverges in D = 5 at six loops, in accordance with 
expectations from previous explicit computations and supersymmetry arguments. This re- 



sult raises various questions about its relation to (2, 0) theory as discussed in ref. [16 
Because (2,0) theory is superconformal, the D = 5 MSYM UV cutoff must be related to the 
gauge coupling, and combining this relation with the S'-duality of the theory compactified 
to D = 4 should lead to strong constraints. Probably the simplest next step is to work out 
the >S-dual extension of the D = 5 two-loop amplitude. 

Through five loops, the planar MSYM four-point integrand has a manifestly nonzero 
behavior in the terms that control the UV divergence in the expected critical dimension D c = 
4 + 6/ L. Thus it is clear, without performing any loop integrals, that the amplitudes diverge 
in the critical dimension. At six loops, we were unable to find an integral representation 
in which all contributions are of the same sign. Therefore we had to explicitly evaluate 
nontrivial integrals in order to answer the question of whether MSYM diverges in D = 5. 
At six loops, practical analytic techniques are not available for generic integrals, so we 
resorted to numerical methods. In particular, we used the sector decomposition method as 
implemented in a modified version of the FIESTA program. 

Our results show that, at least through six loops, the values of the ultraviolet divergences 
in the critical dimension of the planar amplitude approximately follow a simple exponential 
Ansatz. Indeed, extrapolating the results from two, three and four loops using this Ansatz 
matches the calculated values at five and six loops remarkably well. The fact that our 
calculated six-loop value closely matches this extrapolation gives us additional confidence 
that we have computed the six-loop divergence correctly. It also allows us to extrapolate 
the value of the diver gen ces in the critical dimension to even higher loops. Moreover, as 
also noted in refs. [HI, 13[, the ratios of the numerical values of the 1/iV^-suppressed terms 



to the leading-color terms are approximately constant for L = 3,4, 5. Assuming that this 
approximate constancy holds as well for L = 6 gives us a definite prediction for the value 
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of the subleading-color contributions to the divergence in D = 5. While the origin of the 
exponential behavior is still unclear, it does suggest that it might be possible to understand 
the ultraviolet divergences of MSYM in the critical dimension to all loop orders. 

The same integration techniques described in this paper may be helpful for resolving 
other problems. An outstanding question that could be resolved by computation of a higher- 
loop divergence is whether M = 8 and other supergravity theories might be perturbatively 



ultraviolet finite. (For recent reviews see refs. 56J.) The current consensus for M = 8 



supergravity is that a. D = 4 potential counterterm valid under all known symmetries exists 



at seven loops [57J, l58j. (A recent optimistic opinion for all orders finiteness may be found in 
ref. 5{|, while a pessimistic one may be found in ref. (60[.) The same potential counterterm 
could be studied in D = 24/5 at five loops, which should be well within reach of the types 
of integration techniques described here, once the supergravity integrand is constructed. 

Intriguingly, through at least four loops, the explicit values of the Af = 8 supergravity 
divergences (in the same critical dimension as MSYM) are proportional to the same linear 
combination of vacuum integrals as enter the subleading-color divergence of M = 4 theory 
at the same loo p o rder. Moreover, half-maximal supergravity appears to be better behaved 
at three loops [6l[ than had been anticipated [58[. These results emphasize the need for 
further explicit computations at high loop orders, in order to help unravel the ultraviolet 
properties of supergravity theories. Our results demonstrate that evaluations of ultraviolet 
divergences are feasible through at least six loops. Due to the close relation between gravity 
and gauge-theory loop amplitudes [62j], the results presented here also provide a concrete 
initial step towards determining the critical dimension of M = 8 supergravity at six loops. 

As yet, there are no explicit forms of the M = 8 supergravity integrands beyond four loops, 
although recent progress in the nonplanar sector of MSYM at five loops [13J suggests that 
the four-point five-loop amplitude of M = 8 supergravity is within reach. In any case, our 
success at six loops with the sector decomposition method suggests that although difficult, 
an evaluation of the integrals likely to occur at seven loops could be feasible. How hard 
would a numerical evaluation at seven loops be, along the lines discussed here? As we 
discussed, a reasonable guess for the number of sectors of a seven-loop integral is 1.6 x 10 8 . 
If an integral takes 2 seconds, then a 1000-core cluster can evaluate these integrals in about 
3 days. This is perhaps a bit slow as we might have hundreds of graphs and a more complex 
integrand, but with further optimization and a larger cluster even this computation should 
come within reach. It may also be possible to achieve further gains based on converting the 
vacuum i nteg rals to propagator integrals and factorizing them into products of lower-loop 
integrals [2lj], as has been applied recently in maximally supersymmetric theories at four 



and five loops [10|, [13[. There are also other methods for attacking this problem, such as 
the powerful DRA method [63] , which offers much higher precision than can be obtained by 
sector decomposition, provided that an appropriate large system of linear equations can be 
solved symbolically. 

In summary, in this paper we showed that maximally supersymmetric Yang-Mills theory 
diverges at six loops in D = 5, settling the question of whether the link to the (2, 0) theory 
might imply an improved UV behavior. We showed that, even at loop orders as high as six, 
and possibly higher, we can directly determine the UV properties of supersymmetric gauge 
and gravity theories. We also uncovered a simple approximate exponential pattern for the 
values of the divergences in the critical dimension where they first occur. This pattern may 
provide clues toward unraveling the all-loop-order UV structure. It is not obvious how to 
reconcile the appearance of a six-loop divergence in D = 5 MSYM with the finiteness of 
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its UV completion, the (2,0) theory in D = 6. Presumably, the divergence must be cutoff 



by additional degrees of freedom in the UV theory. As discussed in refs. [16|, [18|, there 
are already candidates for these degrees of freedom as nonperturbative states in the D = 5 
theory, so that it may be possible to understand this in D = 5 terms. Perhaps the simplest 
conjecture is that S-dual extensions of the D = 5 amplitudes (as discussed in ref. [l6j]) are 
finite to all orders. 
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